setwd("E:/其他来自老师and同学的数据")
file=read.csv("AllelicEpigenome-sigOnly-AllDocs.csv",header = T)
file$unitID=paste(file$AllelicEpigenome.Variant.properties.Subject.properties.Chromosome.value,file$AllelicEpigenome.Variant.properties.Subject.properties.Position.value,file$AllelicEpigenome.Variant.properties.Subject.properties.Reference.Allele.value,file$AllelicEpigenome.Variant.properties.Subject.properties.Alternative.Allele.value,sep=":")
asm=file

setwd("E:/5hmc_file/2_5hmc_yjp_bam/ASM/")
file=read.csv("20201120/at.least.one.AShM.in.DC.add.BF.beta0.add.CCHC.csv",head=T)
filea=read.csv("20201112做汇总表/all.FDR.sig.at.least.one.add.direction.same.diff.csv",head=T)
filea$id=paste(filea$Chr,filea$Start,sep = ":")
filea1=filea[filea$FDR.sig>1,]

file$id=paste(file$Chr,file$Start,sep=":")
file1=file[file$pattern.not.rm.dupl.num.DC>1,]
file2=file1[file1$BF_in_DC>1,]
file3=file1[file1$BF_in_DC>10,]
> length(intersect(asm$unitID,file3$unitID))
[1] 1
> length(intersect(asm$unitID,file2$unitID))
[1] 11
> length(intersect(asm$unitID,file1$unitID))
[1] 140
> length(intersect(asm$unitID,file$unitID))
[1] 951